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Derivation of effective zero-range one-dimensional (ID) interactions between atoms in tight waveg- 
uides is reviewed, as is the Fermi-Bose mapping method for determination of exact and strongly- 
correlated many-body ground states of ultracold bosonic and fermionic atomic vapors in such waveg- 
uides, including spin degrees of freedom. Odd- wave ID interactions derived from 3D p-wave scatter- 
ing are included as well as the usual even-wave interactions derived from 3D s-wave scattering, with 
emphasis on the role of 3D Feshbach resonances for selectively enhancing s-wave or p-wave scattering 
so as to reach ID confinement-induced resonances of the even and odd-wave interactions. A duality 
between ID fermions and bosons with zero-range interactions suggested by Cheon and Shigehara is 
shown to hold for the effective ID dynamics of a spinor Fermi gas with both even and odd-wave 
interactions and that of a spinor Bose gas with even and odd-wave interactions, with even(odd)- 
wave Bose coupling constants inversely related to odd(even)-wave Fermi coupling constants. Some 
recent applications of Fermi-Bose mapping to determination of many-body ground states of Bose 
gases and of both magnetically trapped, spin-aligned and optically trapped, spin-free Fermi gases 
are described, and a new generalized Fermi-Bose mapping is used to determine the phase diagram 
of ground-state total spin of the spinor Fermi gas as a function of its even and odd-wave coupling 
constants. 
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I. INTRODUCTION 



The spin-statistics theorem, according to which iden- 
tical particles with integer spin are bosons whereas those 
with half-integer spin are fermions, breaks down if the 
particles are confined to one or two dimensions. Real- 
ization of this fact had its origin over forty years ago 
when it was shown [1, 2] that the many-body problem 
of hard-sphere bosons in one dimension can be mapped 
exactly onto that of an ideal Fermi gas, so that many 
properties of such Bose systems are Fermi-like. It is now 
known that this "Fermi-Bose duality" is a very general 
property of identical particles in ID, not restricted to 
the hard-sphere model, and relating strongly interacting 
bosons to weakly-interacting fermions and vice versa. In 
recent years this esoteric subject has become highly rele- 
vant through experiments on ultracold atomic vapors in 
atom waveguides [3-11]. An understanding of their prop- 
erties is important for atom interferometry [12, 13] and 
integrated atom optics [11, 14, 15], which are potentially 
important for development of ultrasensitive detectors of 
accelerations and gravitational anomalies. 

When an ultracold atomic vapor is placed into an atom 
waveguide with sufficiently tight transverse confinement, 
its two-body scattering properties are strongly modified 
and short-range correlations are greatly enhanced. This 
occurs in a regime of low temperatures and densities 
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where both the chemical potential /x and the thermal 
energy fc^T are less than the transverse oscillator level 
spacing hujj_, so transverse oscillator modes are frozen 
and the dynamics is described by a one-dimensional 
(ID) Hamiltonian with zero-range interactions [16, 17]. 
This ID regime has already been reached experimentally 
[8, 18, 19], as has a regime with ji < hwj_ but ksT > Hu)j_ 
[6, 10]. We assume herein that both /j, < frui± and 
ksT < fku± as in [8, 18, 19]. Nevertheless, virtually 
excited transverse modes renormalize the effective ID 
coupling constant via a confinement-induced resonance, 
as first shown for bosons [16, 20] and recently for spin- 
polarized fermionic vapors [21]. At the very low densities 
of ultracold atomic vapors, the 3D interatomic interac- 
tions are usually adequately described by the s-wave scat- 
tering length and a corresponding 3D zero-range pseu- 
dopotential, the input for the original derivation [16] of 
an effective ID interaction between waveguide-confined 
spinlcss bosonic atoms. A more detailed and compre- 
hensive theory has been developed recently [22]. In the 
simplest case of spinless bosons the resultant effective 
ID interaction is of the form gf D 5(z) where g^ D is an ex- 
plicit and nontrivial function [16, 20, 22] of the 3D s-wave 
scattering length, z = z\ — Z2, and z\ and zi are ID coor- 
dinates of the interacting atoms measured along the lon- 
gitudinal axis of the waveguide. This derivation [16, 22] 
will be described in Sec. II A, and in Sec. II B the deriva- 
tion of effective ID pseudopotentials will be presented, 
for application not only to the spinless Bose gas but also 
to the case of spinor Fermi and Bose gases, for which 
the definition of pseudopotentials is much more delicate 
due to wave function discontinuities induced in the zero- 
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range limit by ID odd-wave interactions derived from 3D 
p-wave scattering. In the spatially uniform case (no lon- 
gitudinal trapping potential) the zero-temperature prop- 
erties of an TV-atom spinless Bose gas are determined 
by a dimensionless coupling constant -fs = mgf D /nh 2 
[16, 17, 23] where m is the atomic mass and n = N/L is 
the ID number density, L being the length of the peri- 
odic box. The exact iV-boson ground state was found in 
the spatially uniform case by the Bethe Ansatz method 
in a famous paper of Lieb and Liniger (LL) [23], and 
spawned later development of a powerful and more gen- 
eral approach [24]. An exact solution in the presence of 
longitudinal trapping is not known, but is well approxi- 
mated by a local equilibrium approach [25] . 

Since the density is in the denominator, at sufficiently 
low densities one enters the regime of strong interactions 
and strong short-range correlations where effective-field 
approaches fail, exactly the opposite of the situation in 
3D where the density is in the numerator of the dimen- 
sionless coupling constant. For 75 > 1 the scattering 
reduces to specular reflection (negligible transmission co- 
efficient) and the Hamiltonian reduces to that of impene- 
trable point bosons, whose exact iV-particle ground state 
was found in 1960 by an exact mapping to the ideal Bose 
gas [1], leading to "fermionization" of many properties 
of the Bose gas and related to breakdown of the sym- 
metrization postulate and the spin-statistics theorem [2] . 
This Fermi-Bose mapping method and its application to 
determination of exact N-atom ground states will be de- 
scribed in Sec. Ill A. In Sec. Ill B a very powerful gener- 
alization of this mapping to arbitrary coupling strength 
[26] will be described and its application to determination 
of the ground state of the magnetically trapped, spin- 
aligned Fermi gas will be reviewed, and in Sec. Ill C 
application of two further generalizations of the map- 
ping to determination of the ground state of the optically 
trapped spinor Fermi gas will be described. 
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nji + m 2 being the total mass, yielding H 2 = 
■ Hcom, where 
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where fi — mim 2 /(mi + m 2 ) is the reduced mass, and 
rj_ and Rj_ are the relative and center-of-mass radial co- 
ordinates, respectively The center-of-mass Hamiltonian 
is that of a simple harmonic oscillator whose solution is 
known, hence we focus only on the relative motion of the 
two particles. This reduces the problem to a single par- 
ticle of mass /i, subject to transverse harmonic confine- 
ment, which is scattered by an external potential V(r). 
The central equation which must be solved is therefore 
Schrodinger's equation for the state of relative motion of 
the two atoms 



E-H-V |V(£)> =0 



(4) 



where H = H re i and V is the interatomic potential. De- 
termining the eigenstates of this Hamiltonian, in particu- 
lar within the s-wave scattering approximation for V(r), 
is the central goal of this section. 



2. Definitions and theorems 

In this subsection we briefly review the T-matrix for- 
mulation of scattering theory, which provides a conve- 
nient framework for approaching the present problem. 
Let us first introduce the retarded Green's function for a 
system with the Hamiltonian H and energy E 



II. ATOMIC SCATTERING AND ID 
INTERACTIONS IN TIGHT WAVEGUIDES 

A. Scattering theory in tight waveguides 

1. Formulation of the scattering problem 

We begin from the Hamiltonian for two atoms under 
transverse harmonic confinement and subject to an arbi- 
trary interaction potential 



H 2 



r, 2 



h 2 1 1 

vr2 n vr2 1 2 2 , 1 2 2 

2mi - Vl - 2^ 2 + 2 mi ^ + 2 TO2 ^ 
+ V( ri -r 2 ) (1) 

where m\ and m 2 are the atomic masses, u>± is the trans- 
verse trap frequency, and Vf and rjj_ are the Lapla- 
cian and radial coordinate of the i th atom, respectively. 
This Hamiltonian is separable in relative and center-of- 
mass coordinates R = (miri + m 2 r 2 )/M, r = ri — r 2 , 



Gjj(E) = lim (E + ie-HY 



(5) 



We then define the T-matrix at energy E of the scatter 
V in the presence of the background Hamiltonian H in 
the usual manner as 



[l-Vd 6 (E) 
E \VG 6 (E) 



V 



n=0 



(6) 



the summation form being valid provided that there are 
no difficulties with convergence. 

Two relations on which we will rely heavily are the 
Lippman-Schwingcr relation 



G H+v( E ) 



G^+G^f^^G^E), (7) 



which relates the full Green's function of the system H + 
V to the unperturbed Green's function G^(E) and the 
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T-matrix, and the Lupu-Sax formula [27] 



1-TW25) G 6 (E)-G 6 ,{E) 



H',V 



(8) 
(E), 



which relates the T-matrix of the scatter V in the back- 
ground Hamiltonian H to the T-matrix for the same scat- 
ter but in a different background Hamiltonian H'. 



3. Scattering theory 

In the continuous part of the spectrum of the total 
Hamiltonian H + V its eigenstates can be expressed as a 
sum of an incident and a scattered wave according to 

\1>{E)) = \ME)) + \ME)}, (9) 
where \ipo(E)), the "incident" state vector, satisfies 

Gf(E)\ME))=0, (10) 

we can then express the Schrodinger equation for the 
total system as 



Gf(E) 



V 



(\ME)) + \ME))) = o. (11) 



This equation is readily solved for the scattered wave in 
terms of the unperturbed Green's function and the T- 
matrix, yielding 



\ME))=G 6 (E)f*y(E)\ME)), 



(12) 



which will serve as the basis for our treatment of the 
present scattering problem. 



4- s-wave scattering regime: the reference T-matrix 
approach 

At first glance it may seem that finding the T-matrix 
is no easier than a direct solving of the Schrodinger equa- 
tion (4). We will demonstrate, however, that the T- 
matrix formulation allows for a self-consistent descrip- 
tion of the low-energy part of the spectrum that uses the 
free-space low-energy scattering properties of the inter- 
action potential as the only input. In addition the low- 
energy (s-wave) limit is isolated to a single well-defined 
approximation without requiring the ad-hoc introduc- 
tion of regularization via a pseudo-potential. In this sec- 
tion we first outline this self-consistent low-energy treat- 
ment. We then solve for the T-matrix using the standard 
Huang- Fermi pseudo-potential, showing that the pseudo- 
potential reproduces the exact result in this situation. 

Let the unperturbed Hamiltonian H be a Hamiltonian 
for a single nonrelativistic particle in presence of a trap- 
ping potential U: 



(r|£|V> 



ft 2 V r 2 
2fi 



U(r) 



(13) 



Assume also that the particle is 'perturbed' by a scatterer 
given by 



(r|VW = V(r)(r\ 



(14) 



localized around r = 0. In what follows we will derive 
a low- energy approximation for the T-matrix of the scat- 
terer V in presence of H. It is important to note that, 
by definition, the T-matrix acts only on eigenstates of 
the unperturbed Hamiltonian, which we can safely as- 
sume to be regular everywhere. (This is of course a con- 
straint on the properties of the unperturbed Hamilto- 
nian.) In this case the zero-range s-wave scattering limit 
does not require any regularization of the T-matrix. By 
making use of the Lupu-Sax formula (8), we first de- 
rive the correct form of the T-matrix in the low-energy 
s-wave regime without the introduction of a regularized 
pseudo-potential. In the following section, however, we 
will see that the results we obtain are in agreement with 
the standard Huang-Fermi pseudopotential approach to 
s-wave scattering. 

We begin our derivation by first specifying a "refer- 
ence" background Hamiltonian H' as 



<r|#'|V> = 



h 2 A r 

2/z 



+ E 



(r| 



(15) 



This Hamiltonian is that of a free particle, but with an ex- 
plicit energy dependence included so that the eigenstates 
have zero wavelength at all energies. We note that this 
reference Hamiltonian agrees with the free-space Hamil- 
tonian in the zero-energy limit. While this Hamiltonian 
may seem strange, it is a valid reference Hamiltonian 
which turns out to be useful because the resulting T- 
matrix is energy independent for any scattering poten- 
tial. The Green's function for this Hamiltonian is given 

by 



(r\G 6l (E)\r / ) = - 



1 



27Tft 2 



(16) 



as can be verified by direct substitution into [E — 
H'\ Gfj,(E) = I. In turn the T-matrix of the interaction 

potential V in presence of H' is independent of energy 
and can therefore be expressed as 

(r\f 6l! v(E)\v>)=gD(v,r>), (17) 

where the kernel D is defined as normalized to unity, 



/ 



drdr' D(r,r') = 1. 



(18) 



The normalization coefficient g is then related, through 
the zero-energy scattering amplitude, to the three- 
dimensional scattering length a s according to 



.9 = 



2irh 2 a s 
I 1 



(19) 
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Imagine that the kernel D(r, r') is well localized within 
some radius R. In perturbative expansions at low en- 
ergies this kernel only participates in convolutions with 
slow (as compare to R) functions, in which case it can be 
approximated by a ^-function, 



D(r,r') w <5(r)5(r'). 



(20) 



This straightforward approximation is the key to the s- 
wave scattering approximation. This effectively replaces 
the exact reference T-matrix by its long- wavelength limit, 
so that the reference T-matrix assumes the form 

(r|f A ,^( J B)|r') fe ' fe '« 1/il ^(r)5(rO, (21) 

which is equivalent to 

f A ,p(E)=g\0)(0\, (22) 

where |0) is the position eigenstate corresponding to the 
location of the scatterer. In expression (21) k and k' 
refer to the wavevectors of any matrices which multiply 
the T-matrix from the left and right, respectively. 

If we now substitute the above expression for the ref- 
erence T-matrix into the Lupu-Sax formula (8) for the 
T-matrix under the background Hamiltonian H we ar- 
rive at 

T 6 ,v(E) (23) 

OO 

= ^[.g|0)(0|G^(i?)] ri 3 |0)(0| 

= [l-g(0\G A W\0)+g(0\G A ,W\0)]~\\0)(0\. 

Making use of Eq. (16), we introduce the function x(£)> 
defined as 



X (E) = lim 

r— >0 



<r|^)|0> + ^ 



(24) 



from which wc obtain the following simple expression for 
the T-matrix of the scatterer V in presence of the trap: 



(r|TW(2S)M 



E<^h 2 /fj,R 2 



gS(r) 



1 - 9X(E) 



(r| 



(25) 



From comparing the equations the free-space and 
bound Green's functions obey one can show that the sin- 
gularity in bound Green's function is the same as that in 
the free-space Green's function. Hence, x{E) is the value 
of the regular part of the bound Green's function at the 
origin. 

For the case of transverse harmonic confinement func- 
tion x{E) has been explicitly computed in [22]. It reads 



x(£) = -^C(V2,-( E 



1 



2na± 



))> 



(26) 



where £(s, a) is the generalized Riemann zeta function 
described in the mathematical literature [28] : 



lim 

N^oo 



^ in + a) s 



({s,a) ~ ■ 1 7 - 1 - 1 l - .s (N + a y-- 

Re(s) > 0, -2?r < arg(n + a) < 0. (27) 



Note that no established convention for choosing the 
branch of the irrational power functions exist: the choice 
above is just the most suitable for the needs of this paper. 



5. Effective one-dimensional interaction potential for 
waveguide-confined spinless bosons 

By a direct substitution to the equation for the Green's 
function of the relative motion of two particles in a waveg- 
uide 



(E-H + ie)(v\G 6 (E)\0)=8(v) 



(28) 



it is easy to show that the Green's function can be decom- 
posed to a sum over the transverse modes in the following 
way: 

(t\6 6 (E)\0) (29) 

OO 

= J2( z \Gid(E - hw±{2n + l))|0)^(p)<(0), 

n=0 

where p = xe x + ye y , <p n {p) are the zero- angular- 
momentum eigenstates of the transverse oscillator, and 
Gid{E\d) is the Green's function of a free one- 
dimensional particle: 

{ElD + £S + ie ^ z \ d ^ E - ^ 2n + i))' ) 



= 6(z) 



(30) 



assume now that our energy belongs to the single mode 
window: 



hui± < E < 3hu>±. 



(31) 



In this case all the n^O terms in the expansion (30) ex- 
ponentially decay at large z. Accordingly the solution of 
the scattering problem (9), (12) now resembles a solution 
of a one- dimensional scattering problem: 



ip{r) ' Z = ° {V'o, id(z)+ 



(32) 



(z\Gid(Eid(E))\0)\c/) (0)\ 2 T3d(E)iPo, id (0)} 

where T 3 rj(E) is the strength of the three-dimensional 
T-matrix, i.e., 



t 3D (E) 



gS{r) 



1 - 9X(E) ' 



(33) 



■5 



and Eid(E) = E — fkoj_. Comparing the expression (33) 
and the general form of a scattering solution (9), (12) it 
is natural to interpret the expression in the braces in the 
l.h.s. of (33) as a solution of a one- dimensional scattering 
problem, subject to a scattering potential whose (one- 
dimensional) T-matrix reads 



{z\f 1D (E 1D )\0) - r 1D (E 1D )6(z) 



with 



Tid{E 



ID) 



b (0)\ 2 T 3D (hw ± + E 



ID) 



(34) 



(35) 



For low energies Em <C the strength of the one- 
dimensional T-matrix is approximately 



tid(Eid) 



1 



2a, 



1 



£C(V2) 



■(36) 



The relative error of this approximation scales as 0(k 3 ). 

One can now attempt to introduce an effective one- 
dimensional scatterer whose T-matrix is close or equal 
to the one given above (34), (36). A straightforward 
calculation shows that the T-matrix of a one-dimensional 
delta-potential 



has a form 



V 1D = gf D 5{z) 



n 2 



tid,s{E\d) = 



(37) 



(38) 



where the one-dimensional scattering length af D is re- 
lated to the potential strength by 



9?d = ~ 



(39) 



Comparison of the T-matrices (36) and (38) leads to a 
conclusion that the waveguide scattering T-matrix (36) 
can be exactly reproduced by a delta-potential (37) of a 
scattering length 



,2 r 



2a, 



aj_ 



-C(l/2) 



(40) 



where ((x) is the Riemann zeta function and C(l/2) = 
-1.4603.... 

The above expression reproduces the result of [16] 
where the effective one-dimensional potential has been 
obtained via a straightforward solution of the scattering 
problem. Notice that at a s — \((l/2)\~ 1 a± the coupling 
constant diverges, signifying the so-called "confinement 
induced resonance" (CIR). The significance of this reso- 
nance has been confirmed in ab initio two-body numerical 
calculations with finite-range realistic interatomic poten- 
tials [20] and in many-body Monte-Carlo simulations [29] . 

Further extensions of the zero-range model for the ef- 
fective one-dimensional scatterer can be envisioned. For 



example, if one chooses to to reproduce the scattering 
properties (more precisely the denominator of the one- 
dimensional T-matrix (35)) with a relative error of C(fc 5 ), 
the one-dimensional delta-potential can be replaced by a 
rectangular potential of a finite width 21 and hight /depth 
vq. In the case of repulsive interaction the model poten- 
tial is a rectangular barrier, and for attractive interac- 
tions it is a rectangular well. In the limit of a s <C a± 
their half-widths I are given by 



I 



^ C( 2 3/2) = 0.8081 VoIoT for gf D > 



I = (CO/2)/2) 



2/3 



18 



.06641a, 



for gf D < 



and in both cases the strength vo of the potential is 

vo = \9? D /2l\. (41) 

6. Remarks on p-wave scattering 

Attempts to carry through a program analogous to the 
above for the case of p-wave scattering between polarized 
fermions meet numerous (hopefully technical) obstacles: 
Most of the limiting procedures become mutually nonuni- 
formly convergent and no clear way to identify a correct 
order is visible. In any case the closest candidate for 
the p-wave analog of the free-space three-dimensional T- 
matrix is the pseudopotential introduced in [30] 

(X|T^(£)|V> = ^pW(O) • VV(0) (42) 

that can be shown to reproduce correctly the low-energy 
behavior of the p-wave scattering amplitude. Here V p 
is the p-wave scattering volume, that defines the low- 
energy behavior of the p-wave scattering phase via V p = 
— linifc^o tanS p (k)/k 3 [31]. 

An elegant way around these difficulties has been fond 
recently by Granger and Blume [21], who used a K- 
matrix technique that does not explicitly involve any 
zero-range objects. The analysis of polarized Fermi gases 
presented below is heavily based on the Granger and 
Blume findings. 



B. ID pseudopotentials 

1. Spinless bosons 

This is the simplest case. The ID scattering length 
af D is defined in terms of the ratio of derivative iPb(z) 
and value iPb{z) of the relative wave function just outside 
the range z of the interaction: 

iI>'b( z o) = -</4(-zo) = -{a? D - 2»)"Vb(±zo) , (43) 



which is equivalent in the zero-range limit z 
the familiar LL contact condition [23, 24] 



^(o+) = -<Mo-) 



ft 2 



1>b(0±) 



0+ to 



(44) 
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for the delta function interaction gf D 8(z) provided that 
the scattering length and coupling constant are related 



by gfi 



-ti 2 1 /J,af D where \i is the effective mass to/ 2. 



It follows from the expression for af D derived in [16, 22] 
and Eq.(40) herein that 



9xd = 2a s ftu>_L 



1--|C(1/2)| 



(45) 



implying the existence of a confinement-induced reso- 
nance CIR [16, 20] of the coupling constant as a s is 
tuned via a 3D Feshbach resonance [32] past the reso- 
nance point a s /a± = |C(1/2)| _1 = 0.6848. . .. Hence the 
whole range of ID coupling constants from — oo to +oo 
is experimentally achievable by tuning a s over a narrow 
range in the neighborhood of the ID resonance. It was 
shown recently [20] that this is a ID Feshbach resonance 
between ground and excited transverse vibrational man- 
ifolds. It was shown in [20] and in Sec. II A 5 herein that 
at low longitudinal energies ka± <C 1 the ID scattering 
amplitude generated by the interaction gf D 5(z) repro- 
duces the exact 3D scattering amplitude in the waveguide 
to within a relative error 0(k 3 ). 



2. Spin- aligned fermions 

Consider next a magnetically trapped, spin-aligned 
atomic vapor of spin-i fermionic atoms in a tight waveg- 
uide. The 7V-fermion spin wave function is magnetically 
frozen in the configuration |i • • • |jv ; so the space-spin 
wave function must be spatially antisymmetric, s-wave 
scattering is forbidden, and the leading interaction ef- 
fects at low energies are determined by the 3D p-wave 
scattering amplitude. Such p-wave interactions are usu- 
ally negligible at the low densities of ultracold atomic 
vapors, but they can be greatly enhanced by p-wave Fes- 
hbach resonances [33]. Granger and Blume derived an 
effective one-dimensional K-matrix for the correspond- 
ing two-fermion problem [21] in a tight waveguide. In 
the low-energy [34] domain the K-matrix can be repro- 
duced, with a relative error 0(k 3 ), by the contact condi- 
tion [21, 35] 



V^(0+) - -Vf(O-) = -a( D 1> F (0±) 



(46) 



where 



6V, 



L 1D 



/[i + i2(v4)IC(- 



1/2, l)!]- 1 (47) 



is the odd-wave one-dimensional scattering length, V p = 
a 3 = — lim/j^o t&n 5p(k)/k 3 is the p-wave "scattering 
volume" [31], a p is the p-wave scattering length, and 
C(-l/2,l) = -C(3/2)/4tt = -0.2079... is the Hurwitz 
zeta function evaluated at (—1/2,1) [36]. The expres- 
sion (47) has a resonance at a negative critical value 



Vp c "7 a l = -0.4009 • • •. Note that a b lD -> -oo as V p 
approaches this critical value, implying that the exte- 
rior wave function (i.e., outside the interaction region 



|z| < zq — > 0) satisfies the free-particle Schrodinger equa- 
tion. This is the opposite of the bosonic case, where it 
follows from Eq. (40) that af D — > at resonance. 

In accordance with (46), the low-energy fermionic 
wavefunctions, Eq. (20) of [21], are discontinuous at con- 
tact, but left and right limits of their derivatives coincide. 
V p is tunable via a 3D Feshbach resonance [33] , allowing 
experimental realization of all values of af D from — oo to 
+oo. It will be shown that in this fermionic case the ef- 
fective ID coupling constant is g[ D = — h 2 a^ D /fi 7 which 
can be compared and contrasted with the previously de- 
fined bosonic ID coupling constant gf D = —h 2 /naf D . 
The dimcnsionless fermionic coupling constant is jp — 
mg[ D n/H 2 . Note that the density n is in the numerator, 
whereas it is in the denominator of the bosonic analog 
7b = mg^ D /nh 2 . For 7^ > 1 one has a "fermionic 
TG gas" [35], a fermionic analog of the impenetrable 
Bose gas, called the "Tonks- Girardeau" (TG) gas in re- 
cent literature [14, 18, 25, 29, 37-48]. As previously 
noted, af D — > —00 in the fermionic TG limit, imply- 
ing an interaction-free exterior wave function. It will be 
shown in Sec. Ill B that in this limit the Fermi gas maps 
to an ideal Bose gas, providing a physical explanation of 
the interaction-free nature of the exterior wave function. 
This can be compared and contrasted with the bosonic 
TG gas, which maps to an ideal Fermi gas. 

Although a discontinuity in the derivative is a well- 
known consequence of the zero-range delta function pseu- 
dopotcntial and plays a crucial role in the solution of the 
Lieb-Liniger model [23] , discontinuities of tp itself have re- 
ceived little attention, although they have been discussed 
previously by Cheon and Shigehara [26] and are implicit 
in the recent work of Granger and Blume [21]. For a 
fermionic wave function ipp the discontinuity 2^(0+) 
is a trivial consequence of antisymmetry together with 
the fact that a nonzero odd-wave scattering length can- 
not be obtained in the limit zq — > unless V'f(0±) 7^ 0. 
These discontinuities are rounded off when z > 0, since 
the interior wave function interpolates smoothly between 
the values at z = — zq and z = zq. This is illustrated 
in Fig. 1 for the special case of the fermionic TG gas, 
the limit jf 3> 1 [35]. The potential is chosen to be a 
square well because we will find later that stability of 
the ground state against collapse requires that the corre- 
sponding effective zero-range interaction (zq — ► 0+) be 
negative definite, which can be shown to be the case 
when gf D > and hence both V p < and af D < 0. 
The energy is taken as zero so the exterior solution is 
sgn(z) = ±1; an interior solution fitting smoothly onto 
this is sin(Kz) with k = ^/2/j.Vo/h 2 = n/2z , the crit- 
ical value where the last bound state passes into the 
continuum, a zero-energy resonance. A fermionic con- 
tact condition with a finite scattering length can be ob- 
tained in the limit z — > if k scales with the width Zq 
as k = (n/2z )[l + (2/^) 2 (z /af D )]. 

Following the bosonic case, where the (^-interaction can 
be introduced naturally to cancel the (^-functions result- 
ing from double differentiation of functions with discon- 
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FIG. 1: N — 2 untrappcd fermionic TG gas ground state 
(dashed line) compared with zero-energy scattering solution 
for a square well with range zo and depth Vb corresponding to 
the boundary between no bound state and one bound state, 
a zero energy resonance (solid line), as function of relative 
coordinate z. Units are such that zo = 1. 

tinuous derivatives, in the case of fermions whose wave 
function is discontinuous it is tempting to introduce 5' 
interactions. However, <5'-functions and second deriva- 
tives are ill-defined if used in a convolution with discon- 
tinuous functions. This difficulty is resolved by realiz- 
ing that (as will be shown), the 6' is associated with 
the interior wave function (— z < z < z ) whereas the 
contact condition (46) refers to the z — ► 0+ limit of 
the wave function at z = ±zq 0, just outside the 
range of the interaction potential. Take the Hamilto- 
nian to be H[ D = — (ti 2 /2fi)d 2 + v[ D where v[ D is a 
pseudopotcntial operator to be determined. d 2 is non- 
singular for z 7^ 0, but at the origin there are singular 
contributions. The first derivative is d z tp F {z) — ip F {z ^ 
0) + [-0f(O+) - i> F (0-)]S(z). The second derivative then 
has an additional singular contribution from differentia- 
tion of the delta function: 

d 2 z ^ F (z) = i>" F {z + 0) + [^(0+)-^(0-)]<5'(z) . (48) 

In the second term on the right, S (z) is associated with 
the interior wave function, as we will eventually verify 
explicitly from the z — > 0+ limit of a finite-range inter- 
action, whereas its prefactor arises from the contact con- 
dition just outside the range of the interaction. Define a 
linear regularizing operator d± by 

a±^(z) = (l/2)[^(0+) + ^(0-)] . (49) 

Then the desired pseudopotential is v[ D = g[ D S (z)d±. 
It satisfies a convenient projection property v[ D yb e (z) = 
where ij) e is any even wave function. Although tpp is 
always purely odd (antisymmetric) for the spin-aligned 
Fermi gas, this projection property is useful for the spinor 
(spin-free) Fermi gas to be discussed in Sec. II B 3, where 
the spatial dependence has, in general, both even and 
odd contributions. This is the reason for the general 



form (1/2)[i/j f (0+) + tp F (0— )] of the factor occurring in 
Eq. (49), even though it could be simplified to ^> F (0+) or 
ip F (0—) when acting onip F , for which ip F (0+) = ip F (0—) 
as a consequence of antisymmetry. It is easy to show that 
terms in S (z) cancel from H[ D ip F as a consequence of the 
contact condition (46) provided that g[ D — —h 2 a^ D /fi. 

The physical significance is clarified by starting from 
the same square well already discussed in connection 
with Fig. 1, i.e., v(z) = — Vb when — z < z < z 
and zero when \z\ > zq. The antisymmetric solution 
ipp of the zero-energy scattering equation [(— ft 2 /2/i)<9 2 + 
v(z)]ip F (z) = inside the well is sm(Kz) with k — 
y // 2fiVo/h 2 . The scattering length af D is defined by 

iPf{zq) = -iI>f{-zq) = -(af D - z )ip F {±zo) which is 
satisfied in the limit z — > 0+ if Vb scales with z as 
k = (Tr/2zo)[l+(2/Tr) 2 (zo/a^ D )]. In that limit the bound- 
ary conditions reduce to Eq. (46). Inside the well the ki- 
netic and potential energy terms are —(h 2 /2fi)d 2 ip F (z) — 
— (H 2 K 2 /2iJ,)sin(Kz) and v(z)iPf(z) = — Vb sin(Kz). For 
\z\ < zo, cos(kz) is proportional to a representa- 
tion of 5(z) as z — > 0, since Jf" o cos(nz)f(z)dz — > 

/(°) I-zo cos(nz)dz = f (0)2k" 1 sin(KZo) -» 2z /(0). 
Then its derivative — (k/2zq) sm(nz) is a representa- 
tion of 5 (z). Noting that nz — > n/2 as z — > 
we have — (h 2 /2^i)d 2 ipF(z) = — (fi 2 K 2 /2^) sin(Kz) — > 
(irh 2 /2/i)<5 (z) which agrees with the kinetic energy 
term -(ft 2 /2fj)[^ F {Q+)-ij) F {Q-)}5' (z) since ip F (0+) and 
"0f(O— ) are to be interpreted as ip F (zo) and iI>f{— zq) 
as z — » 0+. Next consider the potential energy 
term inside the well as z — ► 0+: —Vosm(Kz) — > 
— Vb(— 2zo/k)5 (z) — ► (7rfi 2 /2/i)<5 (z). Comparing this 
with vf D ip F (z), using the expression for gf^,, noting that 

Vv(0±) are to be interpreted as V'fC 1 ' 1 ^) for z - * 0, one 
finds that the two expressions for the potential energy 
term agree in that limit. It is clear from this derivation 
that the 5 (z) in the effective ID Hamiltonian can be 
interpreted as the "ghost of the vanished interior wave 
function" , which plays a crucial role in the odd- wave in- 
teraction even in the limit z — ► 0. 



3. Spin-free fermions 

In this section we assume that the spinor Fermi gas is 
optically trapped, so the spins are unconstrained. This 
case is more complicated than the spin-aligned Fermi gas 
even in the absence of explicit spin-spin interactions or 
external spin-dependent potentials, because the require- 
ment of antisymmetry under combined space-spin ex- 
changes (zi,ai) <-> (zj,aj) induces implicit space-spin 
coupling leading to nontrivial spin dependence of the 
wave functions. Here each spin z-component argument o~i 
takes on the values \ or J., or equivalently ±5. Consider 
3D two-body scattering in the waveguide. There are both 
s-wave scattering states, which are space symmetric and 
spin antisymmetric with spin eigenfunctions of singlet 
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form (Tl — IT)i and p-wave scattering states, which 
are space antisymmetric and spin symmetric with spin 
eigenfunctions of triplet form || or J. J. or ^75 (Tl + IT)- 
Assume that the Hamiltonian does not depend on spin. 
Then the spin dependence of fermionic wave functions ipF 
need not be indicated explicitly and they can be written 
as the sum of spatially even and odd parts ?p e and rp . 
The odd part decomposes further into three components 
going with the three triplet spin eigenfunctions, but one 
need not complicate the notation at this point since the 
Hamiltonian acts in the same way on each of them. The 
effective ID interactions are determined by ID scatter- 
ing lengths a\ D for spatially even waves "4> e (z) — ^e{— z) 
related to 3D s-wave scattering and spatially odd waves 
ijj (z) = —ip (—z) related to 3D p-wave scattering. The 
contact condition for ID even-wave scattering is the same 
as the previously-given one (44) with tpB and af D re- 
placed by ip e and a\ D , and the one for ID odd- wave scat- 
tering is Vo(^o) = -ipo{-zo) = -{a°D - Zo)^' {±z ), thc 
zo > version of Eq. (46), with if>p and af D replaced by 
tp and aJ D . In the zero-range limit zq — > 0+ these can 
be combined into [35, 49] 

V/(o+) - v'(o-) = -(a e 1D )- 1 Mo+) + 4>(0-)] 

^(0+) - V(O-) = -a° D W{0+) + V'(0-)] (50) 

where ip(z) = ip e (z) + 4> {z). In this section ip(z) is a 
fermionic function tpp, but we use a more general nota- 
tion because the same equations apply to a spinor Bose 
gas. a\ D and a° D are related to the 3D s-wave scatter- 
ing length a s and the 3D p-wave scattering volume V p by 
Eqs. (40) and (47), with af D replaced by a\ D and af D 
replaced by a° D . Take the Hamiltonian to be 

Hid =-(h 2 /2v)d 2 z +v e 1D +v° 1D . (51) 

Here v° D differs from v[ D of Sec. II B 2 only by the ob- 
vious substitutions, i.e., v° D — g1 D S {z)d± with d± de- 
fined by Eq. (49) with ipp replaced by ip a . The even- 
wave pseudopotential is more complicated than the sim- 
ple delta-function interaction of Sec. II B 1, because the 
delta function is ambiguous at the point z = 0, where 
ip is discontinuous due to the discontinuity in its odd- 
wave component ip a . In order to determine the correct 
form, start from the first derivative d z tp(z) — tp (z 7^ 
0) + [V>(0+) - i>(0-)]S(z) as before. Now the second 
derivative has two singular contributions in addition to 
the nonsingular term ip (z ^ 0), one because in general 
ip (0+) 7^ tjj (0—) and the other from the derivative of 
the delta function: 

d 2 z i>(z) = /(^ ) + [v/(o+)-v'(o-M*) 

+ [^( +)-^(0-)]5'(«) . (52) 

With proper choice of g° D the odd-wave pseudopotential 
v1 D cancels the S (z) term from the kinetic energy, and we 
define the even-wave pseudopotential vf D so as to cancel 



the S(z) term: v\ D — gf D S± where the linear operator 
<5± is defined by 

4V>(2) = (V2M0+) + V(0-Mz) . (53) 

These pseudopotentials satisfy convenient projection 
properties vf D ip — v^ D ^p e = on the even 
and odd parts of V, and their matrix elements are 
(XK D |V>> = ix*(0)[V(0+) + V(0-)] and { X \v° lD W) = 
~hlx (°)]*bP' (°+) + ^'(0-)]. They connect only even 
to even and odd to odd wave functions if we stipulate 
that x(0) = [the average of x(0+) and x(0~ )] if X is 
odd and \ (0) = [the average of \ (0+) and \ (0 — )] if 
X is even. In fact, the wave function and its derivative 
at z = refer to the internal wave function as modi- 
fied by the potential, whereas z = 0+ and z = 0— refer 
to the wave function just outside the range of the po- 
tential, and the above values at z = follow from the 
way the internal wave function interpolates between the 
contact conditions on the exterior wave function. Using 
(50) one finds that terms in 5(z) and S (z) cancel from 
Hid if the even and odd-wave coupling constants are re- 
lated to the scattering lengths by g\ D = —h 2 /na\ D and 

9id = -fr 2 a°i D /V- 

III. FERMI-BOSE MAPPING METHODS AND 
N-ATOM GROUND STATES 

In this section we will review the theory of Fermi-Bosc 
mappings relating the exact TV-particle energy eigen- 
states of systems of fermions and bosons in ID with ef- 
fective zero-range interactions, and application of these 
mappings to determination of the iV-atom ground states. 
This will be done first for the original mapping for im- 
penetrable bosons (7 > 1) [1, 2], then for a very powerful 
generalization to arbitrary values of 7 due to Cheon and 
Shigchara[26], and finishing with a further generalization 
to the case of spinor Fermi and Bose gases, important for 
applications to optically trapped fermionic atoms whose 
spins are unconstrained. 

A. Impenetrable bosons (TG limit) 

It was already pointed out in the famous paper of Lieb 
and Linigcr on the ID Bose gas with delta-function inter- 
actions [23] and in Sees. I and II B 2 that the ID gas of 
impenetrable point bosons is the limit 75 3> 1 of the LL 
gas, the "TG limit" . Tonks gave the first treatment of the 
statistical mechanics of a ID hard-sphere gas [50], which 
was restricted to the classical high-temperature regime 
and provided no information about the extreme quan- 
tum limit characteristic of ultracold atomic vapors. The 
formula for the exact quantum-mechanical ground-state 
energy of the ID hard-sphere Bose gas appeared in a pa- 
per of Bijl where it is quoted without derivation [51], and 
a derivation was published by Nagamiya [52]. Then in 
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1960 one of us [1] and Stachowiak [53] independently red- 
erived the ground-state energy. The Fermi-Bosc mapping 
method was first introduced in [1], although Nagamiya 
had previously noted [52] that in the "fundamental sec- 
tor" z\ < z 2 < ■ ■ ■ < zn the ground state wave func- 
tion of a spatially uniform, ID hard-core Bose gas can 
be written as an ideal Fermi gas determinant, continu- 
ation into other permutation sectors being effected by 
imposing overall Bose symmetry under all permutations 
Zi <-» Zj in spite of the fermionic antisymmetry under 
permutations of orbitals (not coordinates) in the funda- 
mental sector. The mapping theorem is much more gen- 
eral, also holding in the presence of external potentials 
and/or finite two-particle or many-particle interactions 
in addition to the hard core interaction [1]. It also ap- 
plies to the ID time-dependent many-body Schrodingcr 
equation and has been used to treat some time-dependent 
interference properties of the ID hard core Bose gas [54- 
58]. 

We now briefly review the mapping theorem. The 
iV-boson Hamiltonian is assumed to have the structure 
Hid = -(ft 2 /2m) J2jLi ®% + V(z\, - ■ ■ , Zn) where the 
real, symmetric function V contains all external poten- 
tials (e.g., a longitudinal trap potential) as well as any 
finite interaction potentials not including the hard-sphere 
repulsion, which is instead treated as a constraint on al- 
lowed wave functions iPb{zi,- ■ ■ , Zn)'- 

ip B =0 if \zj-z k \<a , l<j<k<N . (54) 

Let ipp{z\, ■ ■ ■ , zn) be a fermionic solution of H\r>i\) = 
Eip which is antisymmetric under all pair exchanges Zj <-> 
Zk, hence all permutations. One can consider ipp to be 
either the wave function of a fictitious system of "spinlcss 
fermions" , or else that of a system of real, spin-aligned 
fermions. Define a "unit antisymmetric function" [1] 

A(z 1 ,---,z N )= Y\ s g n ( z k-zj), (55) 

l<j"<fe<JV 

where sgn(z) is the algebraic sign of the coordinate dif- 
ference z = Zk — Zj, i.e., it is +1(-1) if z > 0(z < 0). For 
given antisymmetric ipp, define a bosonic wave function 
ipB by 

iPb(z\, ■■■,z n )= A(zi, ■■■ , z N )ip F (z 1 , ■■■ ,z N ) (56) 

which defines the Fermi-Bose mapping, ips satisfies the 
hard core constraint (54) if ipp does, is totally symmet- 
ric (bosonic) under permutations, obeys the same bound- 
ary conditions as ipp, and Hip>ipB — EipB follows from 
Hid^Pf = Etpp [1, 2]. In the case of periodic boundary 
conditions (no trap potential, spatially uniform system) 
one must add the proviso that the boundary conditions 
are only preserved under the mapping if N is odd, but 
the case of even N is accomodated by imposing periodic 
boundary conditions on ipp but antiperiodic boundary 
conditions onips- 



The mapping theorem leads to explicit expressions for 
all many-body energy eigenstates and eigenvalues un- 
der the assumption that the only two-particle interac- 
tion is a zero-range hard core repulsion, represented by 
the a — -> limit of the hard-core constraint, the "TG 
gas". Such solutions were obtained in Sec. 3 of [1] 
for periodic boundary conditions and no external poten- 
tial. At the low densities of ultracold atomic vapors it 
is usually sufficient to consider this case, although it has 
been shown by Astrakharchik et al. [29] that for a lon- 
gitudinally trapped LL gas with attractive interaction 
9id < there is a regime within which the equation of 
state is well approximated by taking a > 0. The exact 
ground state is also known in this case [1, 52]. Here 
we limit ourselves to the usual strongly repulsive TG 
limit -fB S> 1- Since wave functions of "spinless" or spin- 
aligned fermions are antisymmetric under coordinate ex- 
changes, their wave functions vanish automatically when- 
ever any Zj — Zk, the constraint has no effect, and the 
corresponding fermionic ground state is the ground state 
of the ideal gas of fermions, a Slater determinant of the 
lowest N single-particle plane- wave orbitals. The exact 
many body ground state was found [1] to have energy 
Eq = (N — N~ 1 )(irhn) 2 /6m where n = N/L is the linear 
number density, and the wave function was found to be 
a pair product of Bijl-Jastrow form 

ipBo = const. | sin^L -1 ^ - Zj)}\, (57) 

where L is the perimeter of the annular trap. In spite 
of the very long range of the individual pair correlation 
factors | sin[7rL _1 (z i — Zj)]\, the pair distribution function 
D(zij), the joint probability density that if one particle 
is found at z% a second will be found at Zj, was found 
to be of short range D(zij) = 1 — [sm(wnzij) / nnzij] 2 . 
Clearly, D(0) = which reflects the hard core nature 
of the two-particle interaction. By examination of the 
excited states the system was found to support prop- 
agation of sound with speed c = irhn/m [1], and it 
was shown that this agrees with the thermodynamic 
formula in terms of the compressibility of the ground 
state. "Fcrmionization" holds only for those properties 
expressible in terms of the configurational probability 
density \ipBo( z i, z 2 , ■ ■ ■ , zn)\ 2 - The momentum distribu- 
tion depends on the single-particle correlation function 
gi(z) = (tp^ (z)ip(O)) (reduced single-particle density ma- 
trix) , which is very different from that of the ideal Fermi 
gas and very difficult to evaluate. Its eigenfunctions 
are plane waves e lkz because of translational invariance 
of the system, and the corresponding eigenvalues define 
the momentum distribution function N(k), the discrete 
Fourier transform of gi(z), the allowed values of k being 
kj = j2ir/ ' L with j = 0, ±1, ±2, ■ ■ -. In a classic tour de 
force Lenard found [59, 60] N(0) to be of order y/N, large 
but much less than the 0(N) value required for Bose- 
Einstein condensation. More generally, gi(z) was found 
[61] to be of order 1/y/k at small k. The corresponding 
momentum distribution is sharply peaked at low k and 
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falls like k~ 4 at large k [37, 39, 62], very different from 
the filled Fermi sea of the ideal Fermi gas, for which N(k) 
is unity for \k\ < &f and zero for |fc| > fcp- 

The exact ground state of the TG gas is also known in 
the presence of a longitudinal trap potential \mui 2 ong z 2 
[63]. It follows from the mapping theorem that the 
exact N-boson ground state ipso is ipBo(zi, • • • , Zn) = 
\ipFo{zi, ■ ■ ■ , zn)\ where tppo is the ground state of N 
spinlcss fcrmions with the same Hamiltonian and im- 
penetrability constraint. The fermionic ground state 
is a Slater determinant of the lowest N single-particle 
eigenfunctions ip n of the harmonic oscillator (HO), where 
ip n (z) = const. e - ^ l 2 H n (Q) with H n {Q) the Hermite 
polynomials and Q = z/z osc , z osc = y/fr/mLJi ong being 
the longitudinal oscillator length. By factoring the Gaus- 
sians out of the determinant and carrying out elementary 
row and column operations, one can cancel all terms in 
each H n except the one of highest degree [64], yielding 
a simple but exact analytical expression of Bijl-Jastrow 
pair product form for the A^-boson ground state: 



V'boCzi, • • -,Zn) = const. 



' N 

it 



-Q-/2 



n \zk-zj\ 

l<j<k<N 

(58) 

with Qi = Zi/z osc . It is interesting to note the strong sim- 
ilarity between this exact ID iV-boson wave function and 
the famous Laughlin variational wave function of the 2D 
ground state for the quantized fractional Hall effect [65], 
as well as the closely-related wave functions for bosons 
with weak repulsive delta-function interactions in a har- 
monic trap in 2D found by Smith and Wilkin [66]. 

Both the single particle density and pair distribu- 
tion function depend only on the absolute square of the 
many-body wave function, and since |^so| 2 = IV'foI 2 
they reduce to standard ideal Fermi gas expressions. 
The single particle density, normalized to N, is n(z) = 

^2 n =o \fn(z)\ 2 an d the P a i r distribution function, nor- 
malized to N(N — 1), is 



D(zi,z 2 ) = n(zi)n(z 2 ) - \A(zi,z 2 ) 



N-l 



A(z 1 ,Z 2 ) = ^ V*n( Z -t)Vn{z 2 ) 



(59) 



n=0 



Although the Hermite polynomials have disappeared 
from the expression (58) for the many-body wave func- 
tion, they reappear upon integrating |?Aso| 2 over (N — 
1) coordinates to get the single particle density n(z) 
and over (N — 2) to get the pair distribution function 
D(z\,z 2 ), and these expressions in terms of the HO or- 
bitals <p n are the most convenient ones for evaluation. 
Some qualitative features of the pair distribution func- 
tion are apparent: In the first place it vanishes at con- 
tact z\ = z 2 , as it must because of impenetrability of the 
particles. Furthermore, the correlation term A (z\,z 2 ) is 
a truncated closure sum and approaches 8{z\ — z 2 ) as 
N — > oo, as is to be expected since the healing length 



in a spatially uniform ID hard core Bose gas varies in- 
versely with particle number [54] . As a result the width 
of the null around the diagonal Qi = Q 2 decreases with 
increasing N, and vanishes in the limit. For \z\ — z 2 \ much 
larger than the healing length, D reduces to the uncorre- 
cted density product n(zi)n(z 2 ), so the spatial extent of 
the pair distribution function is that of the density and 
varies as N 1 / 2 [67]. Detailed gray-scale plots of D(z\, z 2 ) 
in the {Qi,Q 2 ) plane for the cases N — 2, A" = 6, and 
N = 10 are shown in Fig. 1 of [63]. 

The reduced single-particle density matrix with nor- 
malization / pi(z, z)dz — N is 



pi(z,z') = N J ipBo(z 7 z 2 , ■ ■ -,z N ) 

xipBo(z', z 2 , ■ ■ ■ , ZN)dz 2 ■ ■ ■ dzN ■ (60) 

For z ^ z it cannot be expressed in terms of |'0bo| 2 : and 
is therefore very different from that of the ideal Fermi 
gas. The multi-dimensional integral cannot be evalu- 
ated analytically, but in [63] it was evaluated numeri- 
cally by Monte Carlo integration for not too large values 
of N, and grayscale plots are shown in Fig. 2 of [63]. 
More accurate numerical results were found in [68], and 
highly accurate results for large values of N were found 
in [69]. In a macroscopic system, the presence or ab- 
sence of BEC is determined by the behavior of p\(z,z') 
as \z — z'\ — > oo. Off-diagonal long-range order is present 
if the largest eigenvalue of p\ is macroscopic (propor- 
tional to N), in which case the system exhibits BEC and 
the corresponding eigenfunction, the condensate orbital, 
plays the role of an order parameter [70, 71]. Although 
this criterion is not strictly applicable to mesoscopic sys- 
tems, if the largest eigenvalue of p\ is much larger than 
one then it is reasonable to expect that the system will 
exhibit some BEC-likc coherence effects. Thus we exam- 
ine here the spectrum of eigenvalues Xj and associated 
eigenfunctions <j>j{z) ("natural orbitals") of p\. The rel- 
evant eigensystem equation is 



J — ! 



p 1 (z,z')cp j (z')dz' = \j4> 3 {z) 



(61) 



Xj represents the occupation of the orbital <pj, and one 
has Xj — N. Accurate values of the Xj have been 
determined in [69]. In particular, the largest eigenvalue 
Ai was shown to be of order \/~N for large N, as in the 
spatially uniform case. 

Next we examine the momentum distribution, which 
can be shown [63] to be a double Fourier transform of pi : 

/OO /'OO 
dz / dz , p 1 (z,z , )e- lk( - z ' z "> . 
-OO J — oo 

(62) 

The spectral representation of the density matrix then 
leads to N(k) — ^jXj\pj(k)\ 2 where the pj are 
Fourier transforms of the natural orbitals: Pj(k) = 
(27r)- 1 / 2 X^ o 0„(z)e- jfcz dz. The key features are that 
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FIG. 2: Angular cross section versus angle sin(t9) w 6 for 
AT = 10. The dashed line is for the ID gas of impenetrable 
bosons and the solid line is for the corresponding system of 
non-interacting fermions 

the momentum spectrum maintains the sharp peaked 
structure reminiscent of the spatially uniform case [16, 
59, 60] and that the peak becomes sharper with increas- 
ing atom number N. By way of contrast, for a ID Fermi 
gas the corresponding momentum spectrum is a filled 
Fermi sea and can be expressed as N(k) = X^jLi iMj (^)| 2 - 
In a recent paper [72] we devised a scheme to measure 
the momentum spectrum based on Raman outcoupling 
and showned that the angular cross section accurately 
mirrors the momentum distribution. Figure 2 shows the 
angular cross section versus angle for both TV = 10 im- 
penetrable bosons (dashed line) and the corresponding 
system of non- interacting fermions (solid line); see [72] 
for details. 



B. Magnetically trapped, spin-aligned fermions 

Consider first the case N — 2. The spins are frozen in 
the configuration ft by the magnetic field, so the spatial 
relative wave function iPf(z) is antisymmetric (odd in 
z = z\ — z 2 ). The Hamiltonian H[ D and corresponding 
odd-wave pseudopotential were derived in Sec. II B 2. 
Defining a mapped bosonic (even in z) wave function 
by iPb(z) = sgn(z)ipF(z) and mapped scattering length 
af D = af D = aiD where sgn(z) is +1 if z > and —1 if 
z < 0, one finds that ipB satisfies the usual bosonic con- 
tact condition Vs(°+) = "V'bC 0- ) = -( a fr>) _1 ^B(0±), 
the zero-range limit z — > 0+ of Eq. (43). Since the ki- 
netic energy contributions from also agree, one has 
a mapping from the fermionic to bosonic problem which 
preserves energy eigenvalues and dynamics. The relation- 
ship between coupling constants g[ D in H[ D and gf D in 
Hf D = -(h 2 /2^)d 2 z +gf D 6(z) is gf D = h 4 /^g[ D , and by 
(47) this relationship agrees with the low-energy limit of 



Eq. (25) of [21, 34]. In the limit gf D = +oo arising when 
V p —* 0— , this is the N = 2 case of the original mapping 
[1, 2] from hard sphere bosons to an ideal Fermi gas, but 
now generalized to arbitrary coupling constants and used 
in the inverse direction. This generalizes to arbitrary N: 
Antisymmetric fermionic solutions ipp(zi, ■ ■ ■ , zn) are 
mapped to symmetric bosonic solutions Vs(^i) ■ ■ ■ > z n) 
via Eqs. (55) and (56). The Fermi contact condi- 
tions are ipF\ Zj =z t + = -ipF\z 3 =z e - = -{a 1D /2)(d Z] - 
d Ze )ipF\z j =z e± and imply the Bose contact conditions 
(d Zj - d Ze ) : 4'B\z J =z e + = -(d Zj - d Ze )%l) B \z j =z e - = 
— {2/ai]j)ipB\z =zi, and these are the usual LL con- 
tact conditions [23]. The fermionic Hamiltonian is 

H[ D = -(ft 2 /2/x)£jLi^ +9[ D Ei< j <e<NS'(z ji )d ji 

where d jt ip = (l/2)[d Zj ip\ Zj=Ze+ - d Ze ip\ Zj = Ze _}. Al- 
though well-defined in the exact Schrodinger equation 
and in first-order perturbation theory, this fermionic 
pseudopotential becomes ambiguous in higher-order per- 
turbation theory However, after mapping to the bosonic 
Hilbcrt space one has the usual Lieb-Liniger interaction 
gii)S(zji) which is well-behaved in all perturbation orders 
and in second quantization. This generalization of the 
Fermi-Bose mapping theorem, due to Cheon and Shigc- 
hara [26], extends the useful domain of the mapping of 
Eqs. (55) and (56) to the whole range of coupling con- 
stants gf D and g[ D . The first application to the spin- 
aligned Fermi gas is due to Blume and Granger, who 
were led to the mapping by consideration of the zero- 
range limit of a K-matrix formulation [21]. They treated 
only the case N = 2 but did not restrict themselves to 
the low-energy limit considered here. 

The exact ground state [23] of H^ D is known for all 
positive gf D if no external potential or nonzero range 
interactions are present, and the mapping then gener- 
ates the exact TV-body ground state of H[ D - The di- 
mensionlcss bosonic and fermionic coupling constants 
7s = mgf D /nh 2 and -fp — tng^ D n/h 2 introduced in 
Sec. II B 2 satisfy jbJf = 4. The energy per particle 
e is related to a dimensionless function e(-f) available on- 
line [73] via e = (h 2 /2m)n 2 e(j) where 7 is related to 7^ 
herein by 7 = 73 = 4/7^. This is plotted as a function 
of 7_f in Fig. 3. For gir>,F — ► 00 as occurs at a p-wave 
Feshbach resonance, one has the "fermionic TG gas" dis- 
cussed in Sec. II B 2 and [35]. For bosons the TG regime, 
which maps to the ideal Fermi gas, is reached when gf D 
is large enough and/or the density n low enough that 
7b > 1. A similar simplification occurs in the fermionic 
case, where a fermionic TG regime is reached when gf D 
is large enough and/or n high enough that 7f » 1. The 
corresponding fermionic TG gas then maps to the ideal 
Bose gas since 7b7f = 4. 

C. Optically trapped, spin-free fermions 

In this section we assume that the Hamiltonian is spin- 
independent and that the fermionic vapor is trapped in 
a tight optical atom waveguide. The spins are then free 
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FIG. 3: Log-log plot of scaled ground state energy per par- 
ticle e = 2me/h 2 n 2 for the spin-aligned Fermi gas, versus 
dimensionless fermionic coupling constant 7f. 

to assume whatever configuration minimizes the ground- 
state energy. First assume N — 2. The two-body rel- 
ative wave functions are iPf{z) = ip F ( z ) + ip F ( z ) where 
the spatially even part ip e (z) contains an implicit spin- 
odd singlet spin factor, and the spatially odd part ip (z) 
contains implicit spin-even triplet spin factors. States 
of combined space-spin bosonic symmetry can be de- 
fined by the same mapping iPb(z) = sgn(z)tpF(z) used 
in the previous section, which now maps the spatially 
even fermionic function ip F to a spatially odd bosonic 
function ip B and the spatially odd fermionic function ip° F 
to a spatially even bosonic function ip B while leaving 
the spin dependence unchanged. Then the even-wave 
contact conditions for a\ D B follow from the odd-wave 
contact conditions for a° 1D F and the odd-wave contact 
conditions for a° D B follow from the even-wave contact 
conditions for a\ D F . As before, one has a mapping from 
the fermionic to bosonic problem which preserves energy 
eigenvalues and dynamics. The bosonic Hamiltonian is 
of the same form as the fermionic one (51) but with 
mapped coupling constants g\ D B = fr 4 /M 2 ffiD f an< ^ 
9°d b = 1 1 V^9\d f- This generalizes to arbitrary N: 
Fermionic solutions ipF {zi, o\ \ • • • ; zn, <tn) are mapped 
to bosonic solutions i\)b{z\,g\; zm, (Jn) by the usual 
mapping (55), (56), where the spin z-componcnt vari- 
ables <jj take on the values ] and J.. The N-fermion 
and N-boson Hamiltonians are both of the form Hid = 
-(ft 2 /2m)E7 =1 5 z 2 3 +E 

On fermionic states ipF, 9 id an d 9°d are 9 id f arL d 
9 id f> whereas on the mapped bosonic states tpB = Atpp 
they are g e 1D B = h 4 /fi 2 g° DJP and g° D B = h 4 / fi 2 g e 1DF . 
This is discussed in more detail in [49]. 

Assume that both g\ D F > and g° D F > 0. If g^ D F 
is zero then it follows from a theorem of Lieb and Mattis 
[74] that the fermionic ground state has total spin 5 = 
(assuming N even), as shown in the spatially uniform 
case by Yang [75] and with longitudinal trapping by As- 
trakharchik et al. [46] . If g1 D F is not negligible then the 



ground state may not have 5 = 0. In fact, if g\ D F is 
zero then one can apply a theorem of Eiscnbcrg and Lieb 
[76] to the mapped spinor boson Hamiltonian, with the 
conclusion that one of the degenerate Bose ground states 
is totally spin-polarized, has 5 = N/2, and is the product 
of a symmetric spatial wave function ipso and a symmet- 
ric spin wave function. The ground state is then the same 
as the one discussed in the previous section, except that 
now there is an (N + l)-fold directional degeneracy since 
S z can range from —N/2 to N/2. Any 5 = state has a 
higher energy; in fact, for N > 2 the mapped Bose gas is 
partially space-antisymmetric, raising its energy by the 
exclusion principle. 

So far we have considered only the extreme cases where 
either the even-wave or odd-wave coupling constant van- 
ishes. Assume now that they may take on any non- 
negative values. Consider first the case N = 2 of a 
longitudinally trapped spinor Fermi gas, with relative 
spatial wave function iPf(z) and Hamiltonian differing 
from Eq. (51) by addition of a harmonic trap potential 
^fj,ojf ong z 2 . The even and odd- wave coupling constants 
are g\ D F = -h 2 /^ia e 1D F and g° D F = -h 2 a° D F /fi as 
before. ipF may be taken to be either spatially even with 
associated singlet spin function -75 (Tl ~ IT) which has 
5 = 0, or else spatially odd with associated spin function 
which is one of the 5 = 1 triplets ||, J, J., or ^(tl + It)- 
The singlet case has spatially even wave functions iden- 
tical with those of trapped N = 2 bosons. The odd- 
wave pseudopotential then projects to zero so the 5 = 
eigenstates are independent of g° D F , and the even- wave 
pseudopotential reduces to gf D p^( z )- The exact eigen- 
states are known [77], being of the form D„(|£|) where 
D v is a Weber (parabolic cylinder) function [78, 79] and 
z = ^^Jh/2iiuji on g. The absolute value in the argument 
leads to a cusp at z — and the LL cusp condition of 
Eq. (44) and [23] (with g\D,B replaced by g\ D F ) leads 
to a transcendental equation for the allowed values of 
v. T(i — iz/)/T(— if) = —A in terms of the dimen- 
sionless parameter A = g\ D F l2h^/[i/hioi ong . The en- 
ergy eigenvalues are E{y) = (y + ^)huJi on g, the ground 
state is that solution for which v vanishes as A — > 0, 
and its energy Eq is a monotonically increasing func- 
tion of g\ D . Next consider the 5 = 1 (triplet) solu- 
tions, for which ipp is spatially odd. Then the even- 
wave pseudopotential projects to zero, and on carrying 
out the Fermi-Bose mapping tpB(z) = sgn(z)tpF(z) the 
odd- wave pseudopotential is changed to g\ D B 8{z) with 

9id b — h 4 /f J ' 2 9°D f- Thus the 5=1 ground state en- 
ergy is the same function of g\ D B that the 5 = ground 
state energy is of gf D F , and is therefore a monotonically 
decreasing function of g° D F . It follows that the 5 = 
and 5 = 1 ground-state energies are equal on the hyper- 
bola giD, F 9iD,F = h4 /v 2 in the {g\ D F > 9°d,f) plane, 
which forms a phase boundary between the region where 
the absolute ground state has 5 = and that where it 
has 5=1. For g\ D F g° D f < ^/m 2 (below the phase 
boundary) the ground state has 5 = and its energy 



13 



is independent of g° D F , and for g\ D F g1 D F > h / fi 
(above the phase boundary) it has S = 1 and its energy 
is independent of g\ D F . 

Since it depends only on a symmetry argument and 
is independent of LOi ong , it is reasonable to conjecture 
that the above phase boundary is valid for all N and for 
both spatially uniform and longitudinally trapped spinor 
Fermi gases. In other words we make the following 

Conjecture. The ground state of a spinor one- 
dimensional Fermi gas undergoes a quantum phase tran- 
sition at 7 e 7 = 4, being paramagnetic (S = 0) for -f e ^ < 
4 and ferromagetic, with S = N/2 for j e ^ > 4; here j e 
and 7 are the dimensionless even and odd- wave coupling 
constants j e = mgf D F /nh 2 and j Q = mg° D F n/h 2 . 

As a first step in motivating this conjecture, it 
is convenient to change from the representation in 
terms of space-spin antisymmetric wave functions 
iI>f(zi,vi\ ■ ■ ■ ; z n, f jv) to a representation in which states 
with S z = k — T^N are represented by functions 
<Pf(xi, ■ ■ ■ ,x k ;yi, ■ ■ ■ ,VN-k) Ti ■ ■■ Tfeli ■ ■ ■ J-jv-fc, where 
we assume N even and 1 < k < N. tpp is antisymmet- 
ric under permutations of the z-coordinates (x\, ■ ■ ■ , x k ) 
of up-spin atoms and also under permutations of those 
(j/i, • • • , j/jv-fe) of down-spin atoms, but it has no par- 
ticular symmetry under exchanges Xi <-> yj of up-spin 
with down-spin atoms. Formally treating up and down- 
spin particles as different species can be shown [2, 80] 
to be physically equivalent to the usual representation 
iI>f(zi,0i',---',zn,(7n)- The Hamiltonian is similar to 
that of Eq. (3) of [46], but also includes odd- wave in- 
teractions: 

2 / k N-k \ 

\*=1 3 = 1 ) 

V°D,f( x p , x q 

l< P <q<k l<p<q<N-k 
k N-k 

+ Y H [viDA^Vj) + i >°D,F(x il y 3 )} (63) 
»=i j=i 

where the even and odd-wave pseudopotential operators 
are those of Eq.(51) expressed in terms of relative coordi- 
nates x p — x q , etc. Even-wave interactions v\ u F (x p , x q ) 
and vf D F (y p , y q ) between like-spin particles are absent 
because (fF is antisymmetric in both (x\, ■ ■ ■ ,x k ) and 
in (j/i, • • • , yN-k), so the corresponding even- wave inter- 
actions project to zero. However, both even and odd- 
wave interactions between up and down-spin particles are 
present because the states ipp are neither symmetric nor 
antisymmetric under exchanges Xi <-> yj between up and 
down-spin atoms. 

The like-spin odd- wave interactions v° D F (x p , x q ) 
and v° D F (y p , y q ) can be transformed into much 
simpler even-wave interactions by a generaliza- 
tion of the previous Fermi-Bose mapping. Define 
a mapped wave function ips of a two-component 
Bose gas by ^b{xi, ■ ■ ■ , x k ; yi, ■ ■ ■ , Z/jv-fe) = 



A(xi,- ■ ■ , x k )A(y 1 ,- ■ ■ , y N _ k )ip F {x\, ■ ■ ■ ,x k ; y 1 , ■ ■ ■ , y N -k) 
where A is the previously-defined "unit antisymmetric 
function" of Eq. (55), equal to ±1 everywhere, antisym- 
metric in its arguments, and changing sign here only at 
"same-species" collisions x p = x q and y p = y q . Then fs 
satisfies a mapped Schrodinger equation Hb^Pb = Ecps 
with the same eigenvalue E and with 

2 Ik N-k \ 

\i=l j = l ) 

+ 91D,B6(X P ~ X q ) 

l<p<q<k 

+ Y 9iD,B s (y P - y q ) 

l<p<q<N-k 
k N-k 

+ Y^ D A X i >yj)+V° D,F(Xi,y,)] • (64) 
i=l j=l 

Now there are only interactions of LL type [23] be- 
tween like-spin atoms, and these are straightforward to 
treat; they have a mapped coupling constant g\ D B = 
H 4 /n 2 g° D F . On the other hand, the mapping does not af- 
fect the contact condition for collisions Xi = yj of unlike- 
species atoms, so these interactions are the same as in Eq. 
(63) and retain the original fermionic coupling constants 
g\ DF and g° 1D F . 

Consider now the Gross-Pitaevskii (GP) approxima- 
tion, the variational energy with a trial function ip B = 

\Yli=i u t{ x i)]\Ylf=i C u l(yj)}- I n the absence of longitudi- 
nal trapping one has u^{xi) = ui(yj) = L^ 1 / 2 where L is 
the length of the periodic box. When acting on these triv- 
ial wave functions the odd- wave interactions v° D F (xi, yj) 
in (64) project to zero and the even-wave interactions 
vf D F (xi,yj) reduce to simple Lieb-Linigcr interactions 
9iD.F^( x i — yj), so the variational energy per particle in 
the thermodynamic limit (N — > oo,L — > oo, N/L — > n) 
is 

e GP = n[a 2 {gl D , B - gl D , F ) + \{g\ D ^ B + g e 1D , F )] (65) 

where a = S z /N with S z = k- \N. Thus e% p in- 
creases with S z if g\ D B > g\ D F and decreases with S z if 
9 id b < 9 id f- ^ follows that the ground state has S z = 
if 9id,b > 9id,f or equivalent^ g e 1D ^ F g 1D ^ F < h 4 /^ 2 , 
and S = S z = ^ (the maximal value) if g\ D F g\ D F > 
H 4 /fi 2 . This is the same result obtained previously for 
the exact N — 2 ground state in a trap. The ground 
state has S z = if 7 e 7 D < 4 and S = S z = ^ if 7 e 7 > 4, 
so there is a hyperbolic phase boundary 7 e 7 Q = 4 in the 
(7e i 7o) plane. In the S = (a = |) phase the ground 
state energy depends only on gf D B , hence on j Q but not 
on 7e, a result which we will show holds also for the exact 
iV-atom ground state. 

The GP approximation is valid when both gf D F and 
9id b are small enough, more precisely when both j e <C 1 
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and 7 3> 1. Since the product 7 e 7 Q can be made to as- 
sume any desired value without violating these inequal- 
ities, it is reasonable to suppose that the phase bound- 
ary 7 e 7 = 4 holds for the exact ground state. We now 
show that this is true. To see this, note first that since 
ifiB is symmetric in (xi, - ■ ■ ,Xk) and in (t/i, • • • , yjv-fe), 
it follows that when acting on (p B , 9\d b$( x p ~ x q) lii 
equivalent to vf D B (x p ,x q ) defined in connection with 
Eq. (53) withg^ = gf DiB , and similarly gf DiB S(y p -y q ) 
is equivalent to v\ D B (y P , y q ). Furthermore, one may for- 
mally add interactions v° D F (x p , x q ) and v° D F (y p ,y q ) 
since these vanish on tps because of its symmetry in 
(xi, ■ ■ ■ , Xk) and (j/i, • • • , yN~k)- Thus in (64) we may re- 
place g{ D B 5{x p - x q ) by [vl DtB {x p , x g ) + v° D>F {x p , x q )] 

and 9t D ,B S (y P - y q ) b y [«id,b(J/p>2/<z) + ^ DiF (y P ,y q )] 
without changing its action on ip B . After this has been 
done, we note that when g\ D F — g\ D B , the resul- 
tant Hamiltonian reduces to that of a one-component 
Bose gas with both even and odd-wave interactions, with 
ground state energy E (S Z ) which is independent of k, 
hence independent of S z . Hence, on the line 7e7o = 4 
ground states of all values of S z from to N/2 are 
degenerate, and it is easy to show that this remains 
true if spin-independent longitudinal trap potentials are 
added to H B - To complete the proof we need to con- 
vert this into a statement about the total spin S of 
the ground state, not merely its value of S z . To do 
this, first map back to the two-component Fermi gas 

states ip F (xi, - ■ ■ ,x k ;yi, - ■ ■ ,VN-k) Ti ■■■ tfcli ■■■ lN-k- 
They are not in general eigenstates of S, but they can 
be converted into such eigenstates ip F , with the same 
eigenvalues of both energy and S z , by antisymmetriz- 
ing with respect to all combined space-spin exchanges 
{zi.Ui) <-> (zj,crj) after renaming the variables as fol- 
lows: (xi, - ■ ■ ,x k ;yi, - ■ ■ ,yN-k) — > {zi,---,z N ) and (fi 
••• Tfeli ••• lN-k) -» (Ti ••• Tfelfe+i ••• In)- In this 
standard representation the Hamiltonian, total spin, and 
its z-component are mutually commuting, so nondegerate 
energy eigenstates are also eigenstates of S and S z , and 
degenerate ones can be chosen to be such simultaneous 
eigenstates. Let tpFoile, 7o, S z , S) be the lowest such 
state for given values of 7 e and j D . We have shown that if 
7e7o < 4 then such a ground state has S z — 0, if 7 e 7o > 4 
it has S z = ^r, and if 7 e 7 Q = 4 it can have any value of S z 
from to -j- If 7e7o > 4 this state also has S = ^ since 
S z = y implies S = y. If 7 e 7o < 4 then the ground 
state also has S = because if it had S > than there 
would exist states with < \S Z \ < S of the same energy, 
contradicting the monotonic dependence of energy on S z 
demonstrated in (65), which we assume to hold also for 
the exact ground state. It follows that the hyperbola 
7e7o = 4 is the boundary between a ground state S = 
phase (7e7o < 4) and aS = | phase (7 e 7o > 4). 



IV. PROSPECT AND CRITIQUE 



It follows from the above that the dependence of the 
ground-state energy on total spin S is more complicated 
than envisioned in either [46] or [49]: It is indeed true 
that if 7 is exactly zero, then the ground state has 5 = 
as assumed in [46], but if 7 Q ^ then by approaching 
the confinement-induced resonance (CIR) of ■y e implied 
by Eq. (45), one can in principle always make 7e7o > 4 
even if < 7„ <S 1, thereby inducing a phase transition 
from the paramagnetic S = phase to the ferromagnetic 
S = N/2 phase. In the opposite case j S> 1 encountered 
near the CIR of j Q implied by Eq. (47) a phase transition 
in the opposite direction (S = N/2 to S = 0) would in 
principle occur if 7 e — > 0+, but in practice there is no 
way of making j e that small. In the region 7 e 7 Q > 4 
where S = N/2 the exact ground state is totally spatially 
antisymmetric and hence spin-aligned, so its energy is 
independent of -y e and given by the results in [49] and 
Sec. II B 2 herein. In the region 7 e 7o < 4 where S = the 
ground state is thus far only known for the case j = 0, 
where it was determined by Yang [75] in the spatially 
uniform case and by Astrakharchik et al. [46] in the 
longitudinally trapped case. No analytical or numerical 
results are yet known for the ground-state energy in the 
region 7 e 7 Q < 4 if j D ^ 0, but it should be investigated by 
numerical calculations. These will be more complicated 
than the previous ones [46] due to the presence of both 
even and odd-wave interactions between up and down- 
spin atoms. However, if the ground state is real and 
nodeless [81] they should be feasible, perhaps by using an 
interaction potential consisting of a narrow and deep well 
to represent the odd- wave interaction (see Sec. II B 2) 
together with a somewhat broader "soft rod" potential 
to represent the even-wave interaction. 



Our proof of the exact ground-state phase boundary 
7e7o — 4 is a "physicist's proof" and we make no claim 
of mathematical rigor. In the first place, as pointed out 
in Sec. II A 6, a rigorous derivation of the zero-range, 
ID limit of the odd- wave interaction between fermions in 
tight waveguides does not exist at present, although we 
believe that Eq. (47) is a correct zero-range, low-energy 
consequence of the K-matrix treatment of Granger and 
Blumc [21]. Furthermore, for N > 2 our proof of the 
phase boundary is not completely rigorous since we had 
to assume that the ground state energy is a single- valued 
function of S z both for 7 e 7 < 4 and for 7 e 7 > 4. We 
feel that this is justified since we proved it to be true 
in the GP approximation and showed that the criteria 
for validity of the GP approximation can be satisfied on 
the phase boundary 7 e 7 D = 4. However, a truly rigorous 
proof of this phase boundary does not yet exist, and we 
hope that one will be forthcoming. 



15 



Acknowledgments 

We are very grateful to Doerte Blume for helpful com- 
ments and for communications regarding her closely- 
related works with Brian Granger [21] and with As- 
trakharchik et al. [46], and to Ewan Wright for refer- 



ences [50-53]. We thank Michael Moore and Thomas 
Bergeman for their invaluable assistance in preparation 
of the Les Houches lecture notes [22] on which Sec. II is 
based. This work was supported by Office of Naval Re- 
search grant N00014-03-1-0427 (M.D.G. and M.O.) and 
by NSF grants PHY-0301052 and PHY-0070333 (M.O.). 



[1] M. Girardeau, J. Math. Phys. 1, 516 (1960); for the case 
of nonzero hard core diameter see footnote 7 (a). 

[2] M.D. Girardeau, Phys. Rev. 139, B500 (1965), particu- 
larly Sees. 2, 3, and 6. 

[3] M. Key et al, Phys. Rev. Lett. 84, 1371 (2000). 

[4] D. Miiller et al, Phys. Rev. Lett. 83, 5194 (1999). 

[5] J.H. Thywissen, R.M. Westervelt, and M. Prentiss, Phys. 
Rev. Lett. 83, 3762 (1999). 

[6] K. Bongs et al., Phys. Rev. A 63, 031602 (2001). 

[7] J. Denschlag, D. Cassettari, and J. Schmiedmayer, Phys. 
Rev. Lett. 82, 2014 (1999). 

[8] M. Greiner et al, Phys. Rev. Lett. 87, 160405 (2001) and 
Nature 415, 39 (2002). 

[9] J. Arlt and K. Dholakia, Phys. Rev. A 63, 063602 (2001). 
[10] A. Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001) and 

Nature 415, 39 (2002). 
[11] N.H. Dekker et al., Phys. Rev. Lett. 84, 1124 (2000). 
[12] P.R. Berman, ed., Atom Inter ferometry (Academic Press, 
Boston, 1997). 

[13] R. Dumke et al, Phys. Rev. Lett. 89, 220402 (2002). 

[14] J. Reichel and J.H. Thywissen, arXiv:cond-mat/0310330 
(2003) (to be published in Proceedings of Euro Sum- 
mer School on Quantum Gases in Low Dimensions, Les 
Houches, April 2003). 

[15] J. Schmiedmayer, Eur. Phys. J. D 4, 57 (1998). 

[16] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[17] D.S. Petrov, GV. Shlyapnikov, and J.T.M. Walraven, 
Phys. Rev. Lett. 85, 3745 (2000). 

[18] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 
Rev. Lett. 91, 250402 (2003). 

[19] B.L. Tolra et al, arXiv:cond-mat/0312003 (2003). 

[20] T. Bergeman, M. Moore, and M. Olshanii, Phys. Rev. 
Lett. 91, 163201 (2003). 

[21] B.E. Granger and D. Blume, arXiv:cond-mat/0307358 
(2003). 

[22] M.G Moore, T. Bergman, and M. Olshanii, arXiv:cond- 
mat/0402149 (2004) (to be published in Proceedings of 
Euro Summer School on Quantum Gases in Low Dimen- 
sions, Les Houches, April 2003). 

[23] E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 

[24] V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Func- 
tions (Cambridge University Press, Cambridge, 1993). 

[25] V. Dunjko, V. Lorent, and M. Olshanii, Phys. Rev. Lett. 
86, 5413 (2001). 

[26] T. Cheon and T. Shigehara, Phys. Lett. A 243, 111 
(1998) and Phys. Rev. Lett. 82, 2536 (1999). 

[27] Adam Lupu-Sax, Quantum Scattering Theory and Appli- 
cations, PhD Thesis, Harvard University (1998). 

[28] T. M. Apostol, Introduction to Analytic Number Theory 
(Springer, 1976). 

[29] G.E. Astrakharchik, D. Blume, S. Giorgini, and B.E. 
Granger, arXiv:cond-mat/0308585 (2003). 



[30] R. Roth and H. Feldmeier, Phys. Rev. A 64, 043603 

(2001) . 

[31] H. Suno, B.D. Esry, and C.H. Greene, Phys. Rev. Lett. 

90, 053202 (2003). 
[32] J.L. Roberts, NR. Claussen, S.L. Cornish, EA. Donley, 

EA. Cornell, and GE. Wieman, Phys. Rev. Lett. 86, 

4211 (2001). 

[33] GA. Regal, C. Ticknor, J.L. Bonn, and D.S. Jin, Phys. 

Rev. Lett. 90, 053201 (2003). 
[34] This follows from the low-energy limit of Eq. (19) of [21] 

which includes transverse zero-point energy hco± , so zero 

longitudinal energy implies E — hu)±. 
[35] M.D. Girardeau and M. Olshanii, arXivxond- 

mat/0309396 (2003). 
[36] ET. Whittaker and G.N. Watson, A Course of Modern 

Analysis (Cambridge University Press, 1952), pp. 

265-269, particularly the top equation on p. 269, which 

applies to the case here where the first argument 

of the Hurwitz zeta function is negative. See also 

www.itu.dk/bibliotek/encyclopedia/math/h/h393.htm 

and references therein. 
[37] M. Olshanii and V. Dunjko, Phys. Rev. Lett. 91, 090401 

(2003). 

[38] P. Ohberg and L. Santos, Phys. Rev. Lett. 89, 240402 

(2002) . 

[39] M. Olshanii and V. Dunjko, New Journal of Physics 5, 
98 (2003). 

[40] D.M. Gangardt and GV. Shlyapnikov, Phys. Rev. Lett. 
90, 010401 (2003) and New Journal of Physics 5, 79 

(2003) . 

[41] G.E. Astrakharchik and S. Giorgini, arXivxond- 

mat/0212512 (2002). 
[42] N.M. Bogoliubov, C. Malyshev, R.K. Bullough, and J. 

Timonen, arXiv:cond-mat/0306735 (2003). 
[43] E.H. Lieb, R. Seiringer, and J. Yngvason, Phys. Rev. 

Lett. 91, 150401 (2003). 
[44] K.V. Kheruntsyan,D.M. Gangardt, P.D. Drummond, and 

GV. Shlyapnikov, Phys. Rev. Lett. 91, 040403 (2003). 

arXiv:cond- mat/0308219 (2003). 
[45] P.D. Drummond, P. Deuar, and K.V. Kheruntsyan, 

arXiv:cond-mat/0308219 (2003). 
[46] G.E. Astrakharchik, D. Blume, S. Giorgini, and L.P. 

Pitaevskii, arXiv:cond- mat/0312538 (2003). 
[47] F. Gerbier, arXiv:cond- mat/0401507 (2004). 
[48] GP. Berman, F. Borgonovi, F.M. Izrailev, and A. 

Smerzi, Phys. Rev. Lett. 92, 030404 (2004). 
[49] M.D. Girardeau and M. Olshanii, arXivxond- 

mat/0401402 (2004). 
[50] L. Tonks, Phys. Rev. 50, 955 (1936). 
[51] A. Bijl, Physica 4, 329 (1937), Note II. 
[52] T. Nagamiya, Proc. Phys. Math. Soc. Japan 22, 705 

(1940). 

[53] H. Stachowiak, Acta Univ. Wratislaviensis 12, 93 (1960). 



16 



[54] M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 84, 
5691 (2000). 

[55] M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 84, 
5239 (2000). 

[56] M.D. Girardeau, K.K. Das, and E.M. Wright, Phys. Rev. 

A 66, 023604 (2002). 
[57] K.K. Das, G.J. Lapeyre, and E.M. Wright, Phys. Rev. A 

65, 063603 (2002). 
[58] K.K. Das, M.D. Girardeau, and E.M. Wright, Phys. Rev. 

Lett. 89, 170404 (2002). 
[59] A. Lenard, J. Math. Phys. 5, 930 (1964). 
[60] A. Lenard, J. Math. Phys. 7, 1268 (1966). 
[61] H.G. Vaidya and C.A. Tracy, Phys. Rev. Lett. 42, 3 

(1979). See also corrected expansions in M. Jimbo, T. 

Miwa, Y. Mori, and M. Sato, Physica D 1, 80 (1980). 
[62] A. Minguzzi, P. Vignolo, and M. Tosi, Phys. Lett. A 294, 

222 (2002). 

[63] M.D. Girardeau, E.M. Wright, and J.M. Triscari, Phys. 

Rev. A 63, 033601 (2001). 
[64] A.C. Aitken, Determinants and Matrices (Oliver and 

Boyd, Edinburgh and London, 1951), p. 112. 
[65] R.B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
[66] R.A. Smith and N.K. Wilkin, Phys. Rev. A 62, 061602 

(2000). 

[67] E. B. Kolomeisky,T.J. Newman, Joseph P. Straley, and 
Xiaoya Qi, Phys. Rev. Lett. 85, 1146 (2000). 

[68] G.J. Lapeyre, Jr., M.D. Girardeau, and E.M. Wright, 
Phys. Rev. A 66, 023606 (2002). 

[69] P.J. Forrester, N.E. Frankel, T.M. Garoni, and N.S. 
Witte, Phys. Rev. A 67, 043607 (2003). 

[70] O. Penrose and L. Onsager, Phys. Rev. 104, 576 (1956). 



[71] C.N. Yang, Rev. Mod. Phys. 34, 694 (1962). 
[72] M.D. Girardeau and E.M. Wright, Phys. Rev. Lett. 87, 
050403 (2001). 

[73] e(7) is made available online by V. Dunjko and M. Ol- 
shanii at http://physics.usc.edu/~olshanii/DIST/. 

[74] Elliott Lieb and Daniel Mattis, Phys. Rev. 125, 164 
(1962). 

[75] C.N. Yang, Phys. Rev. Lett. 19, 1312 (1967). 
[76] Eli Eisenberg and Elliott H. Lieb, Phys. Rev. Lett. 89, 
220403 (2002). 

[77] M.A. Cirone, K. Goral, K. Rzazewski, and M. Wilkens, 
J. Phys. B 34, 4571 (2001). Note, however, that their 
Fig. 3 is incorrect, since the curves for c / do not have 
the necessary cusps at z — ; for example, the c = oo 
curve (TG limit) should have its left half reflected up, 
yielding an even function with the necessary cusp. 

[78] P.M. Morse and H. Feshbach, Methods of Theoretical 
Physics, Vol. II (McGraw-Hill, New York, 1953), pp. 
1565-1567 and 1641. 

[79] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Se- 
ries, and Products (Academic Press, New York, 1980), 
pp. 1064-1067. 

[80] A. Messiah, Quantum Mechanics (North-Holland Pub- 
lishing Co., Amsterdam, 1962), Vol. II, pp. 619 ff. 

[81] This has been proved by Lieb and Mattis [74] assuming 
that the interactions are ordinary potentials. However, 
our odd-wave interaction involves both a derivative and 
another limiting procedure, and we are not aware of a 
proof applying to this case. 



